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A COMPUTATIONAL EXAMINATION OF 

DIRECTIONAL STABILITY FOR SMOOTH AND CHINED FOREBODIES AT HIGH-ct 

R. Ravi 

William H Mason 

Virginia Polytechnic Institute and State University 
SUMMARY 

Computational Fluid Dynamics (CFD) has been used to study aircraft forebody flowfields at 
low-speed, high angle-of-attack conditions with sideslip. The purpose is to define forebody 
geometries which provide good directional stability characteristics under these conditions. The 
flows over the experimentally investigated F-5A forebody and Erickson forebody, previously 
computed by the authors, were recomputed with better grid topology and resolution. The results 
were obtained using a modified version of cfl3d to solve either the Euler equations or the 
Reynolds equations employing the Baldwin-Lomax turbulence model with the Degani-Schiff 
modification to account for massive crossflow separation. Based on the results, it is concluded that 
current CFD methods can be used to investigate the aerodynamic characteristics of forebodies to 
achieve desirable high angle-of-attack characteristics. An analytically defined generic forebody 
model is described, and a parametric study of various forebody shapes was then conducted to 
determine which shapes promote a positive contribution to directional stability at high angle-of- 
attack. An unconventional approach for presenting the results is used to illustrate how the positive 
contribution arises. Based on the results of this initial parametric study, some guidelines for 
aerodynamic design to promote positive directional stability are presented. 
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INTRODUCTION 


Current and future fighter aircraft are likely to be required to operate at high angles-of-attack 
where they experience flowfields that are dominated by large regions of separated vortical flows. 
Considerable research is being done both in the experimental and computational areas to 
understand the physics of such complex flows. A good understanding of these flows would enable 
engineers to design fighter aircraft to achieve better maneuverability at high angles-of-attack. At 
high angles-of-attack the forebody aerodynamic characteristics make significant contributions to the 
complete configuration aerodynamics. The surveys by Chambers (Ref. 1) and Chambers and 
Grafton (Ref. 2) present the basis of the current understanding of high angle-of-attack 
aerodynamics. 

One of the specific aerodynamic characteristics of interest is directional stability. For the F- 
5A, which has good high angle-of-attack characteristics, it has been shown experimentally (Ref. 3) 
that the forebody makes a significant positive contribution to directional stability at angles-of-attack 
above which the vertical tail ceases to be effective. That forebody had a smooth cross section, 
although it was not axisymmetric. The current authors recently demonstrated that the experimental 
results could also be predicted using Computational Fluid Dynamics (CFD) methods (Refs. 4 and 
5). The ability to reproduce previously obtained experimental results meant that it would be valid to 
use CFD to try to design shapes for specific aerodynamic characteristics at high angle-of-attack, 
where large regions of separated flow are present. 

Future advanced fighters are likely to possess chine type forebodies, as evidenced by the YF- 
22 and YF-23 configurations. For these aircraft, high levels of agility are demanded, and the 
aerodynamic characteristics at high angle-of-attack play an important role in determining aircraft 
handling qualities and agility. 

Because of the interest in chine-shaped forebodies, a key issue in the application of 
computational methods to forebody design is the ability to treat chine sectional shapes. Few general 
chine-shaped forebody wind tunnel tests are available to use for comparison with computational 
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methods. One is the wind tunnel investigation conducted by Erickson and Brandon (Ref. 6, the 
“Erickson Forebody”). In that study the chine effects were investigated for a generic fighter 
configuration, and pressure distributions were measured on the chine forebody. All forebody 
results were acquired in the presence of the wing. More recently Kegelman and Roos (Refs. 7 
and 8) studied experimentally the influence of cross sectional shape on the vortex flowfield at high 
alpha. They compared the surface pressures and the aerodynamic loads between a circular, 
elliptical and a chined cross section at high angles of attack. Hall (Ref. 9) studied the influence of 
the forebody cross sectional shape on wing vortex-burst location. This study also involved the 
comparison of a two chine cross sections with a circular section. 

The results were obtained using cfl3d (Ref. 10) to solve either the Euler equations or the 
Reynolds equations employing the Baldwin-Lomax turbulence model with the Degani-Schiff 
modification to account for massive crossflow separation. Version 1.1 of the code with the 
modifications as described in Ref. 5 was used in all the computations. 

In this report we first repeat the results obtained on the F-5A forebody (Ref. 5) using a grid 
strategy better suited to the geometry to assess possible grid sensitivity on the previous results. 
Secondly, we compare computed predictions with the experimental data for the chined class of 
forebodies using the “Erickson Forebody” (Ref. 6) at a = 30° (5° and 10° sideslip)and at a = 40° 
(10 sideslip). The above two cases were used to establish confidence in the methodology being 
used for the analysis of various generic cross sectional forebodies. 

An analytic model which can be used to systematically study forebody aerodynamics for 
families of forebody shapes at high alpha is proposed. Using this model, a computational study is 
earned out to determine which shapes lead to the best directional stability characteristics. The 
reference parameters used in computing the forces and moments for the cases studied in this report 
are presented in Table 1. The report concludes with some guidelines for high angle-of-attack 
forebody design. 
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wall friction velocity, / p 
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inner law variable, yu /v 
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P angle of sideslip, deg 

v kinematic viscosity 

p density 

x w shear stress at the wall 

ACp difference between leeward and windward Cp across the vertical plane of 
symmetry 
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F-5A FOREBODY 


The wind tunnel experiment demonstrating the dominant contribution of the F-5A forebody 
to directional stability at high angle-of-attack was simulated computationally in the first part of this 
work (Ref. 5). This forebody had been tested by Sue Grafton, et.al. at NASA Langley Research 
Center and the results are available in Ref. 3. The geometry math model and the comparison with 
the wind tunnel model was described in detail in Ref. 5. In that study the grid was constructed 
from two dimensional O-type cross flow grids which are longitudinally stacked, constituting a 
single block H-0 topology as shown in fig. 1. It is difficult to resolve the flow details near the 
nose using an H-O topology. Hence, we investigated the same geometry using an alternate grid 
system to assess possible grid effects on the results. 

F-5A Grid Details 


The inviscid calculations on the F-5A (Ref. 5) were repeated on the new grid shown in 
fig. 2. Ibis grid consists of two blocks, where the first block used a C-O topology to improve the 
grid resolution at the nose. This grid was generated using a transfinite interpolation grid generator 
provided by Ghaffari (Ref. 11). The first block extends from the nose to the point where the flat 
sidewall starts i.e., 14.025 inches from the nose, as explained in Ref. 5. The inviscid calculations 
were performed on a grid which used 32 axial, 93 circumferential and 45 radial points ( 32 x 93 x 
45 ). The outer boundary extends 32.7 inches radially outward and is comparable to the length of 
the forebody which was 31.02 inches. The second block used the previous H-0 grid topology 
with 13 axial, 93 circumferential and 45 radial points (13 x 93 x 45). The C-O grid generator used 
for the first block requires a user specified normal distance to the first grid point and the distance of 
the outer boundary as the input. The H-0 grid generator used for the second block uses the 
distance of outer boundary and a stretching parameter as the input. Care was taken to ensure that 
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the distance of the first grid normal to the surface is the same for both the blocks at the interface. 

Figures 3 and 4 show the grid used for inviscid calculations at different cross sections 
downstream from the nose. Figure 3(a) shows the entire cross sectional grid at FS 14.02 and fig. 
3(b) shows the details near the body at the same station. Figures 4(a) and 4(b) contain the same 
information at FS 29.61. Because of the presence of the flat sidewall at sections downstream, the 
grid points were clustered near the maximum half breadth points forward of the flat sidewall. This 
provided adequate definition of the flat wall portion of the forebody. 

A grid refinement study was done for both inviscid and turbulent solutions for an angle-of- 
attack of 40°. The grids used in this study were the same H-O grids used in the first part of this 
work (Ref. 5). The baseline inviscid grid had 33 (axial), 93 (circumferential) and 45 points in the 
radial direction. The baseline viscous grid had 33 (axial), 93 (circumferential) and 65 points in the 
radial direction. During the grid refinement study, the number of points in the radial direction were 
increased with improved radial stretching, so that at least four fine grid points were present in the 
first grid point of the crude grid. The circumferential and axial densities were kept the same. The 
inviscid refined grid had 90 points in the radial direction while the refined viscous grid had 100 
points radially. 

Results and Discussion of Computations on the F-5A Forebody 


Inviscid calculations were performed for a = 30° and ft = 5° to compare the results of this 
new grid system with those obtained using a H-O grid earlier in Ref. 5. The boundary conditions 
were the same in both the cases except on the axis that runs from the nose to the upstream farfield 
boundary where a singularity type boundary condition was imposed for the new grid. In the earlier 
computations this boundary was a part of the surface and so an inviscid boundary condition was 
imposed. 
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Figures 5 and 6 show the comparison of the inviscid surface pressures between the two grid 
systems at FS 6.58 and FS 26.77 respectively. It is very difficult to identify a difference between 
the two results at these stations. For the two-block grid system, FS 6.58 is in the first block where 
the surface grid was different from that of the H-O grid used earlier in Ref. 5. FS 26.77 is in the 
second block, where the surface grid was same as that used earlier. A comparison of the 
longitudinal variation of C p on the leeward plane is shown in fig. 7. This figure shows the 
improved resolution of the solution near the nose with the new C-0 grid. There was negligible 
change in the value of the directional stability C n p After demonstrating that the two grid systems 
produced similar results, a grid resolution study, as well as a detailed study of the solutions, was 

done on the single block H-0 grid used in Ref. 5. 

Figure 8 shows the sign convention used in computing the side force and yawing moment. 
The F-5A forebody experimental directional stability data from Ref. 3 are shown along with the 
computed inviscid and viscous results in fig. 9. The computed results revealed the same trend 
found in the wind tunnel data and were already presented in Refs. 4 and 5. Additionally, we show 
the results obtained with the refined grid for both inviscid and turbulent cases at a = 40°. 
Although the results changed slightly with grid resolution, the trends were the same in both the 
cases. 

Figure 10 shows the axial distribution of side force contributing to the yawing moment 
presented in fig. 9 at a = 40°. The importance of the viscosity in producing the positive stability is 
clearly demonstrated. The viscous solution develops a significant restoring force, with a positive 
side force over most of the forebody and generally increasing with downstream distance. This is a 
consequence of the increasing asymmetry of the forebody vortices with distance from the nose. 
The inviscid solution shows essentially no side force over the majority of the forebody. Figures 1 1 
and 12 provide the circumferential pressure distributions at two stations for both inviscid and 
turbulent cases. The corresponding cross sectional shape, the direction of incoming flow and the 
origin of reference for the angular measure are shown below each of these figures. The negative 
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peak pressures are due to the vortices on the upper surface of the cross section and are shown more 
clearly in the following flow visualization pictures to be presented in fig. 14. The asymmetry in the 
pressure distribution due to the sideslip can be seen in fig. 11, and is much more noticeable in 
fig. 12. At FS 14.02 the viscous solution results clearly show the effect of the vortices, with two 
low pressure regions, denoted B and C, underneath the vortices. The low pressure peaks A and D 
are due to the attached flow accelerating around the highly curved sides of the body. At station FS 
29.61 the inviscid results contain four distinct low pressure peaks corresponding to the high 
curvature regions at the cross section comers. Considering viscous effects, the turbulent flow is 
separated and the primary vortices are moving away from the body, as shown later in fig. 14. The 
small low pressure peak at C in fig. 12 is due to the primary windward vortex. The primary 
leeward vortex is sufficiently far off from the surface and hence the suction created by it is 
insignificant 

Figure 13 contains the pressure differences, AC p , between the leeward and windward sides 
of the body at the same stations at which the pressures were plotted in figs. 1 1 and 12. These 
provide insight into the distribution of side force at a particular station to help explain the effect of 
viscosity in creating the restoring force. Although the viscous effects are primarily associated with 
the vortex and separated flowfield on the top side of the forebody, the effects of viscosity are seen 
to alter the balance of pressures between the sides of the body over most of the side projection. It is 
particularly interesting to notice that the near zero side force associated with the inviscid flow arises 
as a delicate balance between a side force in one direction on the lower portion of the body, and a 
side force in the opposite direction on the upper part. The effects of viscosity are to reduce the 
magnitudes of the vortex suction peak levels as well as producing a shift in the location of the 
ACp curve which results in a distribution which has a much larger net side force. 

Figures 14(a) and 14(b) show the cross sectional stagnation pressure contours at axial 
stations x = 14.02 inches and x = 29.61 inches from the nose for the viscous calculation at the 
same flow conditions shown in figs. 11 and 12. The incoming flow is the same as shown in 
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figs. 11(b) and 12(b). The leeward (LHS) vortex is farther away from the surface than the 
windwaid (RHS) vortex. The asymmetry of the low pressures on the body under the vortices is 
generally considered to be pulling the body to smaller sideslip, and thus provides a stabilizing 
moment However, we have shown in fig. 13 that the side force is affected by the separated flow 
indirectly through its effect on the pressure distribution over virtually the entire surface. These local 
effects of the vortices actually act primarily on the essentially flat top-surface and don’t directly 
contribute to the side force. At FS 29.61 we can also see the secondary vortices under the primary 
vortices. 

Figures 15 and 16 show the inviscid and turbulent calculation pressures at FS 14.02 and FS 
29.61 plotted as vectors perpendicular to the surface. In these diagrams, the surface is treated as a 
line of zero pressure and the vectors going outward from the surface represent negative pressure 
coefficients. These diagrams should be studied in conjunction with pressure plots of fig. 1 1 and 
12. At FS 14.02 the viscous diagrams show clearly the effect of the vortices resulting in two low 
pressure peaks on the upper surface. At FS 29.61, the inviscid results show two peaks as the flow 
accelerates around the sides to the upper surface. The viscous results shown in fig. 16(b) can be 
explained by looking at the flow structure in fig. 14(b). The leeward vortex is sufficiently far off 
from the surface and hence the leeward vortex suction peak is almost insignificant On the 
windward side, as expected, the windward vortex is closer to the surface and the low pressure 
peak is therefore still visible. 

Figure 17 shows the vortex path development along the body. The leeward vortex (here on 
the RHS of the body) is seen to be rising above the body much faster than the windward vortex 
(LHS). The windward vortex is actually “blown back” over the forebody, and is moving along the 
top surface near the center. This is also evident from fig. 14. 
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ERICKSON CHINE FOREBODY 


Erickson and Brandon experimentally investigated the chine effects on a generic fighter 
configuration and have published the detailed pressure data over a large range of angles of attack 
and sideslip (Ref. 6). This flowfield on such a model was computationally investigated in the 
earlier part of this work and the details were presented in Ref. 5. The geometry math model and the 
comparison with the wind tunnel model are also described in detail there. In that study, the grid 
was constructed from two dimensional O-type cross flow grids which are longitudinally stacked, 
constituting a single block H-0 topology as was done earlier in the case of F-5A. Here we 
investigate the same geometry using an alternate grid system and also study the grid resolution 
requirements for a chined forebody. 

Erickson Forebody Grid Details 


The inviscid calculations on the Erickson forebody were repeated on the alternate C-O grid 
shown in fig. 18. The baseline inviscid calculation grid had 45 points in the radial direction and 
101 points in the full circumferential direction. Longitudinally, the grid was clustered near the nose 
with 25 stations on the forebody as shown in fig. 18. The axial grid planes were defined at 
stations corresponding to the experimental measured stations. These were at a distance of 7.19, 
13.56 and 19.94 inches from the nose along the length of the body. The smoothing of the surface 
unit normals introduced some grid skewness near the chine nose as well as around the chine edge. 
This was done to avoid large cell volume discontinuities. 

As compared to the inviscid solution grid, the viscous calculation used a grid with 65 points 
in the radial direction, and with longitudinal and circumferential grid points remaining identical 
with the grid used for the inviscid calculations. The baseline grid was established with sufficient 
normal clustering near the surface to adequately resolve the laminar sublayer in the turbulent 
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boundary layer flow. This grid produced an average normal cell size of approximately 10 -4 /. At 
the ffeestream conditions used in the computations for the Enckson forebody (Af,*, — 0.2, Re/ — 
1.02 x 10 6 based on model length, and a = 20° ) the baseline grid typically resulted in a value of 

y + « 2 at the first mesh point above the surface. 

Figure 19 shows the grid used for inviscid calculations at the last section downstream from 
the nose. Figure 19(a) shows the entire cross sectional grid at FS 14.02 and figs. 19(b) and 19(c) 
provide the details near the surface and chine edge respectively. 

A grid refinement study was done with both the inviscid and turbulent grids. In each of these 
cases the n um ber of grid points were doubled in the normal direction with increased clustering in 
the normal direction. The circumferential and axial densities were kept the same. Approximately 
four fine grid points were packed in the first cell of the baseline grid for both the fine inviscid and 
the fine turbulent grids. The fine Navier-Stokes grid provided ay + value of approximately 0.5. 

Results and Discussion of Computations on the Erickson Forebody 

Inviscid calculations were performed for a = 30° and p = 0° to compare the results of this 
new grid system with those obtained earlier using an H-0 grid in Ref. 5. The H-0 grid used earlier 
had 33 axial stations with 25 on the surface and 8 ahead of the nose. With the new C-0 topology 
grid strategy, the radial and circumferential grid densities were kept the same. As in the case of 
F-5A, the boundary condition on the axis that runs from the nose to the upstream farfield boundary 
was altered to a singularity type boundary condition. In the earlier computations this boundary was 
a part of the surface and so an inviscid boundary condition was imposed. 

Figures 20 and 21 show the comparison of the inviscid surface pressures between the two 
grid systems at FS 7.19 and FS 13.56 respectively. The surface grid for the C-0 grid topology 
was changed to accommodate more points near the chine edge while keeping the circumferential 

*The experimental Mach number of .08 was not used because of previous problems using other 
Navier-Stokes codes. An examination of the results showed that the largest Mach number in the 
flowfield was well removed from the compressibility range. 
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grid density same as that of the H-O grid. The difference is almost insignificant as was seen in the 
case of F-5A. The advantage of using this grid system is that you can maintain the same grid 
density on the surface while eliminating the upstream grid extension. 

Figures 22 - 24 present the computed upper surface pressure distributions at three stations 
obtained on the isolated forebody along with experimental data on the forebody-wing model for 
various angles of attack and sideslip. The details of the experimental investigation are available in 
Ref. 6. Some of the results arising from earlier computations are reproduced from Ref. 5 in these 
figures. The surface grid used in later computations were stretched to accommodate more points 
near the chine edge. Figure 22 shows the upper surface pressures for the a = 30° and 5° 
case. At the section closest to the nose (FS 7.19) the inviscid computations predict the pressures 
very close to the experimental values. At stations further downstream the agreement deteriorates. 
At FS 19.94 the wind tunnel data appears to reflect the higher local incidence induced by the wing 
flowfield. The inviscid refined grid results show a suction peak in better agreement than the 
baseline grid at the first station, but produce little change further downstream. Turbulent viscous 
effects do not change the pressure levels at the mid section of the forebody, but do have some 
effect on the peak suction pressure level. The peak suction pressures were reduced, as expected, 
resulting in poorer agreement with the experimental data. In the turbulent flow case the refined grid 
solution resulted in only minor changes in the pressure distribution. The trend remains the same 
when the sideslip is increased to 10° as shown in fig. 23. The fine grid Navier-Stokes calculation 
was not done for this case. Figure 24 shows the pressures for a = 40° and /J = 10°. Erickson 
and Brandon ( Ref. 6 ) suggest that the extent of upstream influence of vortex breakdown 
occurring downstream was found to differ at different combinations of angles of attack and 
sideslip. For example, vortex core bursting had occurred on the windward side at a = 40° 
whenever the sideslip angle exceeded 5°. The computations do not include a model of the wing 
effects. We did not observe any vortex bursting in the solution. 
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Figure 25 shows the side force computed for the Erickson forebody. Both the inviscid and 
the viscous solutions show similar trends, and the minor grid effects indicate that the solutions are 
grid resolved. Here, in contrast to the smooth forebody cross section results for the F-5A, both 
inviscid as well as the turbulent results develop restoring forces, with a positive side force over 
most of the forebody and generally increasing with downstream distance. This is expected because 
of the fixed separation lines along the edge of the chine, regardless of viscosity, and is in marked 
contrast to the smooth cross section results obtained on the F-5A forebody (Ref. 4). There the 
inviscid and viscous solutions were completely different, with the inviscid solution providing 
essentially no side force. The vortical flow in this case is being governed essentially by inviscid 
phenomena. The directional stability characteristics in fig. 26 show the stabilizing effect of the 
chined forebody over the entire range from 20° to 40° angle-of-attack. Qualitatively, the trend 
shown by both Euler and Navier-Stokes grids are very similar. This observation is important, and 
provides a basis for deciding on the solution strategy to be used for the parametric computations on 
a generic forebody to be discussed later.The directional stability computed for this forebody is 
similar using either Euler or Navier-Stokes solutions at 30° angle of attack. At a = 40°, the refined 
Navier-Stokes grid calculation resulted in improved correlation with Euler results. 

Figures 27 and 28 show the inviscid and turbulent calculation pressures at FS 7.19 and FS 
19.94 plotted as vectors perpendicular to the surface. As before, the surface is treated as a line of 
zero pressure and the vectors going outward from the surface are proportional to the negative 
pressures. These diagrams should be studied in conjunction with pressure plots of fig. 23. Unlike 
the case of F-5A, the inviscid and turbulent cases are very similar at both the stations because of 
the fixed separation line as discussed earlier. The flow decelerates as it approaches the chine edge 
because of the change of the body cross sectional shape, and separates at the chine edge. 

Figures 29 and 30 show inviscid and turbulent stagnation pressure contours respectively at 
two different stations. These clearly show the chine-edge generated vortices. The position and 
magnitude of the primary vortices are nearly identical in both the inviscid and turbulent cases. The 
turbulent solution also shows the formation of secondary vortices near the chine edge due to 
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boundary layer separation. For chine shapes the effects of viscosity are a secondary effect on 
vortex size, position and strength. Strong vortex formation can be seen all along the forebody in 
fig. 31 with the leeward vortices rising above the surface much faster than the windward vortices. 
Such strong vortex formation on bodies with sharp chines is responsible for positive directional 
stability even at 20° angle of attack which was not found in the F5-A case. 


SOLUTION STRATEGY FOR PARAMETRIC FOREBODY GEOMETRY STUDY 


Based on the analysis of the computational solutions obtained on the Erickson chine 
forebody, a solution strategy for forebody shaping study was chosen. When /J was fixed at 5° it 
was shown in the case of the Erickson chine forebody that the inviscid pressures were very close 
to the experimental data and the side force and Cnp trends were qualitatively similar and nearly the 
same for the Euler and turbulent flow computations. Though refining the grid made a slight 
improvement in the Euler results, it was very expensive considering the minor change in the 
results. Hence, it was decided that to assess aerodynamic trends arising from forebody geometry 
variations on chine-shaped forebodies, the computations could be done with an Euler analysis and 
the baseline grid. 

To study the advantage of using multigrid and multisequencing, the inviscid flow over a 
generic analytical forebody was computed at a = 30° and /J = 5°. Three levels of sequencing 


were used with multigridding on each level. The surface pressures as shown in fig. 32 were 
identical when the residual went down to the same order of magnitude in both cases. However, 
there was a 33% reduction in CPU time. After this approach was established, the remaining Euler 
calculations were performed with three levels of sequencing and multigridding on each level of 
sequencing. 
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ANALYTIC CHINE FOREBODY STUDY 


To study geometric shaping effects on forebody aerodynamic characteristics, an analytical 
forebody model with the ability to produce a wide variation of shapes of interest was defined in 
Ref. 5. This generic forebody model makes use of the equation of a super-ellipse to obtain the 
cross sectional geometry. The super ellipse, used previously to control flow expansion around 
wing leading edges (Ref. 12), can recover a circular cross section, produce elliptical cross sections 
and can also produce chined-shaped forebodies. Thus it can be used to define a variety of different 
cross sectional shapes. 

The super-ellipse equation for the forebody cross section was defined in Ref. 5 as: 

(if- 

Kb) 

where n and m are adjustable coefficients that control the surface slopes at the top and bottom 
plane of symmetry and chine leading edge. The constants o, and b correspond to the maximum 
half-breadth and upper or lower centerlines respectively. Depending on the value of n and m, the 
equation can be made to meet all the requirements specified above. The case n = m = 0 
corresponds to the standard ellipse. The body is circular when a = b. 

When n = -1 the sidewall is linear at the maximum half breadth line, forming a distinct crease 
line. When n < -1 the body cross section takes on the cusped or chine-like shape. The derivative 
of ilb with respect to y / a is: 



(2 + ni 

\ 

dz _ 

U + n 

) 

dy 

[ 1 — H 

r l+iA 

<2+n) 


where I = z lb and y = y/ a. Asy 1, the slope becomes: 
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00 n >- 1 

0 n <-l 

-(2 + m)y l+m n = —l 

Different cross sections can be used above and below the maximum half-breadth line. Even 
more generality can be provided by allowing n and m to be functions of the axial distance x, 
although in this study the parameters n and m were taken to be constants with respect to x. The 
parameters a and b are functions of the planform shape and can be varied to study planform 

effects. Notice that when n = -1 the value of m can be used to control the slope of the sidewall at 
the crease line. 

Using the generic forebody parametric model defined above, and the computational strategy 
developed based on the Erickson forebody results, an investigation of directional stability 
characteristics of various chine-shaped forebody geometries was made. It was decided to analyse 
the effect of changing bia , chine angle and combinations thereof. This range of cross sectional 
shapes provides an extremely broad design space to investigate aerodynamic tailoring of forebody 
characteristics through geometric design. 

For the present study the following cases were initially selected : 

(a) Geometrical parameters: 

m = 0 

-1.5 <n< -1.0, An = -0.25 
0.5 < b/a < 1.5, Ah/a = 0.5 

(b) Flow conditions: 

20° £ a < 40°, A a= 10° 

0°<fi< 5°, A/J= 5° 

The resulting cross sectional shapes are shown in fig. 33. The computational study was 
earned out to determine the shape which leads to the largest increase in directional stability. This 


dz__ I 
dy~\ 
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test case matrix, shown in Table 2, resulted in 54 different configurations with symmetrical upper 
and lower surfaces, showing how large the possible set of cases could be without careful selection. 
The p = 0° cases were initially included to compare the flow physics with and without sideslip. 
However, with C n = 0 at 0 = 0° and the number of cases being excessive, the p = 0° cases were 
eliminated. Further combinations were eliminated as the study progressed and the results 
examined. Some asymmetric upper/lower cross section geometries were also analysed. These 
geometries were defined using different b/a or different n for upper and lower surfaces. 

It was also decided that the planform shape would initially be defined to be similar to the 
Erickson chine case and to study the effects of varying cross section geometry. In this calculation 
the moment center for the computation of the directional stability was kept fixed at the value used in 
the Erickson forebody test (Table 1). Based on the best cross sectional shape, limited planform 
effects were studied. The total CPU time used for the Euler study is given in Table 3. 

Discussion of Results for the Generic Chine Forebodies 
Effect of varying b/a 

This study was conducted for cross sectional shapes with m = 0andn = -1.5 and bla = 0.5, 
1.0 and 1.5 (see fig. 33d). Figure 34 shows C n ^v s angle of attack with b/a as the varying 
parameter. It is interesting to note that the contribution to positive directional stability increases as 
b/a decreases at a fixed angle-of-attack. b/a = 0.5 is the best cross section in promoting positive 
directional stability. An understanding of these results requires an examination of the flowfield 
details presented below. 

Figure 35 shows the variation of the side force with the axial distance at each angle-of-attack. 
Near the nose the force is initially destabilizing, being negative for all cases computed. Moving aft 
from the immediate vicinity of the nose, the trend is reversed and the side force starts to increase 
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toward positive values. The side force becomes more positive with increasing angle-of-attack. In 
general, the side force becomes increasingly negative as the value of b/a increases, making the 
body more unstable. However, some crossover occurs at the aft end of the body at the higher a, 
where the bla = 0.5 case is not as positive as the b/a = 1 case. 

Figures 36 to 38 show the A Cp vs z plots at a typical cross section ( x = 18.35 ). The 
integration of this pressure difference produces the side force values presented in the fig. 35. The 
cross section below the chine edge always makes a negative contribution to the side force. Above 
the chine edge there is an abrupt large positive spike in the side force. This arises because of the 
asymmetry in strength and position of the vortices. At a= 20° the shallow b/a = .5 case produces 
a much larger spike than the b/a = 1.5 case. At higher a the b/a = 1 case has nearly the same size 
spike. 

The asymmetry in the position and strength of the windward and leeward vortices which is 
responsible for the positive side force on the forebody is shown in fig. 39 for a = 30° and b/a = 
0.5 and 1.5. Figure 39(a) shows the minimum static pressure found in the vortex over the length 
of the body. In this case the lower pressure for the b/a = 0.5 geometry is much stronger compared 
to the b/a =1.5 case. Also, the windward vortex for this geometry is much stronger than the 
leeward vortex resulting in a larger asymmetry. This corresponds to the large difference in 
directional stability shown in fig. 34. In the sideview shown in fig. 39(b), for b/a = 0.5 both the 
vortices are further away from the chine line than in the b/a = 1.5 case, and they are above the top 
centerline, allowing communication between the windward and leeward vortices. In the planform 
view, fig. 39(c), the b/a = 0.5 case shows more lateral movement particularly in the aft region than 
the b/a = 1.5 case. Here the windward and leeward vortices are separated by the large hump on the 
upper surface all along the length of the forebody, and thus restricts the influence of one vortex on 
the other, as weU as the vortex movement. This is illustrated in fig. 40, which presents stagnation 
pressure contours to show the increase in vortex movement as b/a decreases 

Using these results the physics of chined forebody aerodynamics emerges. A shallow upper 
surface (b/a = 0.5) results in a stronger, more asymmetric vortex system compared to a deep 
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surface {bla = 1.5). A deep lower surface results in a larger negative contribution to directional 
stability. Hence, higher bla for the upper or the lower surface is undesirable. 

Effect of varying chine angle 

In this study bla was held constant at 0.5 (corresponding to the best result obtained above) 
and n was varied over -1.5.-1.25 and -1.0, which increases the edge angle from a sharp chine to a 
sharp edge ( Fig. 33a ). Recall that theoretically the chine edge has a zero angle when n = -1.5 and 
n = -1.25 and therefore has a 180° slope discontinuity. When n = -1.0 the included edge angle is 

finite (127°) and the slope discontinuity is smaller. 

The effect of changing the shape parameter n on the directional stability is shown in fig. 41. 
Essentially, all the results are similar at a = 20° and 30° but show differences at or = 40 . The 
sudden decrease in C n ^ for n = -1.0 at a = 40° was further investigated by looking at the side 
force variation in fig. 42. Based on the results shown in this figure for the n = -1 case over the 
axial distance from about 3 to 23, the source of the decrease of at a = 40° forn = -1.0 can 
be identified. This result provides an indication of how to keep C n ^ from becoming too positive 
at high angles-of-attack. Figures 43 to 45 show the &C p vs z plots at a typical cross section 
( x = 18.35 ). At a = 20°and 30° the effect of the chine angle is predominant on the upper surface. 
Though the behavior changes on the upper surface, the area under the curves remains nearly the 
same. At a = 40° the area under the curve suddenly decreases for the n — -1.0 case and this leads 
to a decrease in side force at this cross section. Figure 46 shows the vortex strength and position 
for the case of n = -1.0. This shows that the side force could arise from the asymmetry in both the 
relative strengths and relative positions of the windward and leeward vortices. 

The vorticity being generated due to flow separation has been shown to be proportional to the 
square of the velocity at the separation point in Ref. 14. When n < -1 the slope discontinuity is 
maximum at the chine edge, and results in large velocities approaching the separation point. This 
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results in larger vorticity being generated at the chine edge for these cases. Figure 47 shows the 
square of velocity at the separation point plotted for different chine angles at a = 40°. The 
n = -1.0 case is distinctly different than the other cases. When n < -1 the edge angle is zero and 
hence the strengths of the corresponding leeward and windward vortices are comparable. Also, 
very close to the nose the leeward vortex is stronger than the windward vortex leading to a negative 
side force. As the axial distance increases the vorticity shed on the windward side increases and 
hence the side force is positive. Such observations were also made by Kegelman and Roos based 
on experimental results in Ref. 7. When n = -1, as expected, the vorticity shed is much less and of 
an entirely different character because of reduced slope discontinuity. Moving downstream from 
the nose, the edge with the largest separation velocity switches sides several times. This is reflected 
in the side force plot of fig. 42(c). In this case, very close to the nose the windward vorticity shed 
is larger than leeward vorticity leading to a positive side force. As we move aft, the side force 
changes sign as the relative shed vorticity strength changes. 


Effect of un symmetrical bla 


Unsymmetrical cross sections were generated using different values of b!a for the upper and 
lower surfaces while keeping the same functional form with m = 0 and n = -1.5. This maintains the 
zero chine edge angle for all the cases. Two cases were tested. The first one had bla = 0.5 for top 
and bla =1.5 for bottom. The second one had bla = 1.5 for top and bla = 0.5 for bottom. Figure 
48 shows the cross sectional shapes together with the computed C n ^ for these bodies alongside 
the results already presented for symmetrical bla. The Erickson forebody result is also included, 
which is geometrically similar with symmetrical bla lying between 0.5 and 1.0. The shallow upper 
surface is seen to provide higher C n ^ than the shallow lower surface geometry. This is because 
the shallow upper surface results in a stronger vortex and provides a bigger contribution to stability 
than the use of a shallow lower surface to reduce the negative contribution to stability. 
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Using the b/a = .5 case as the baseline in fig. 48, it is interesting to contrast the effects of 
increasing the body height above and below the chine line. The rate of C n ^ reduction at a fixed 
angle of attack due to additional body height above the chine line is less than the rate of 
reduction when the height is added below the chine line. In addition, the change of C n ^ with 
angle of attack differs depending on whether the height is added above or below the chine 
line.Here too, the loss in rate of change with angle of attack is less when thickness is added above 
the chine line instead of below it. 


Effect of unsymmetrical cross sections to vary chine angle 

Unsymmetrical cross sections were generated using different values of shape parameter n for 
the upper and lower surfaces while keeping the same bla = 0.5 which was found to be the best 
ratio earlier. Such a variation of n would vary the chine angle. The effect of varying this parameter 
on the directional stability is shown in fig. 49. The chine angles were zero for symmetrical cross 
sections with n <- 1 and were finite for all other cases shown in that figure. Only the symmetrical 
case with n = -1.0 which had the highest chine angle shows a sudden decrease in C n ^ at a = 

40°. This difference in behavior with the different chine angles suggests the existence of a critical 
angle which controls the rate of feeding of the vortex as the angle-of-attack changes. 


Effect of varying the planform shape 


The planform shape for the forebodies studied thus far was same as that of the Erickson 
forebody. This planform is shown in fig. 50. The parameter XN shown for the tangent ogive 
forebodies is the distance from the tip of the nose to the station where the planform span becomes a 
constant. The side force variation in figs. 35 and 42 showed that most of the positive side force 
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came from the aft portion of the forebody where the chine line was swept nearly 90°. Hence it was 
postulated that expanding to a constant cross section faster would give greater positive side force. 
Because the Erickson planform approximates a tangent ogive with XN = 18, the alternative 
planform was chosen to expand faster with XN =7, as shown in fig. 50. 

The effect of the planform variation on the directional stability is shown in fig. 5 1 . There is a 
small increase in C„ for a fixed cross sectional shape with bla = 0.5, m = 0 and n = -1.0. One 
other cross section, with a flat lower surface, was computed with this planform, and resulted in a 
C np incr ease. This supported our previous assertion that a smaller bla on the lower surface 
reduces the adverse contribution to C n at a = 20° and also at a = 40°. Here, note that the chine 
included angle is much less than the symmetrical case. The directional stability continues to 
increase at a = 40°, rather than remain nearly constant, reinforcing the idea that a critical chine 
angle might exist which reduces extreme contributions to stability at high angle-of-attack. 

Figures 52 and 53 show the effect of planform shape on side force variation at a = 20° and 
ot — 40 respectively. As expected, after the initial negative side force, the rate of increase of side 
force is greater in the aft portion of the forebody for the blunt nosed planform. Also note that at 
a = 40 , the double hump is eliminated with a blunt-nosed planform and with a flat bottom 
surface the configuration is even better. However very close to the nose the side force is more 
negative. A look at the slopes and curvatures of the different planforms in fig. 54 shows that the 
tangent ogive planform has a large negative curvature close to the tip of the forebody. 

CONCLUSIONS AND DESIGN GUIDELINES 


A number of conclusions arise based on the results obtained here. For chined-shaped 
forebodies, where the separation position is not influenced by viscosity, the Euler solutions were 
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found to be in reasonably good agreement with the results of Navier-Stokes calculations using the 
Baldwin-Lomax turbulence model as modified by Degani and Schiff. Thus Euler solutions could 
be used to carry out the parametric study. CFD has been used to explicitly identify the method in 
which the pressure distribution on the chine contributes to the directional stability. An 
unconventional approach to presentation and evaluation of forebody aerodynamics has been 

introduced. 

For aerodynamic design consideration the following guidelines were obtained: 

• The best ratio of maximum half-breadth to the maximum centerline width proves to be 
b/a = 0.5 among the cases analysed for positive directional stability. In general, lower b/a 
for both the upper and lower surfaces improves directional stability. In cases where higher 
b/a is a requirement, it is better to increase the lower surface b/a which results in a smaller 
penalty than if we were to increase upper surface b/a. The rate of change of C np * s a ^ so a 
function of the surface to which the thickness is added. 

• The effect of chine angle on the directional stability characteristics was found to be 
insignificant except when the chine angle was large. There could be a critical chine angle 
beyond which it becomes an important factor (we did not attempt to find one in this study). 
If such a critical angle exists, it provides an indication of how to keep ^ from becoming 

too positive at high angles of attack. 

• The positive contribution to the stability is seen to come from the aft portion of the 
forebody where the chine line is swept nearly 90°. Changing the planform shape by 
allowing it to expand faster to a constant value increases the C n ^ only by a small amount. 
However, the behavior of the side force plots vary significantly for different planform 
shapes. 
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Reference Parameters 

F-5A 

Erickson 

Generic 

Mean Aerodynamic Chord c 

16.08 in 

32.04 in 

32.04 in 

Wing Span b’ 

52.68 in 

46.80 in 

46.80 in 

Model Length / 

31.025 in 

30.00 in 

30.00 in 

Reynolds Number Re [ 

1.25 x 10 6 

1.02 x 10 6 

1.02 x 10 6 

Reference Area S re jr 

754.56 in 2 

1264.32 in 2 

1264.32 in 

Moment Reference Center from Nose 

57.72 in 

12.816 in 

12.816 in 


Table 1. Reference Data Used in Computing Forces and Moments 
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Matrix of Cases for “Symmetric” Chine Forebodv Directional Stability' 
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Table 2. Total Cases for Parametric Study 


Each inviscid “crude grid” run = 3400 CPU seconds + 200 sec = 3600 sec 


effect 


planforms hl£L s 


«’s £s 


n s 


total 


bla 1 

sideslip (/3 = 10°) 1 

chine angle (extra) 1 

split bla 1 

split chine angles 1 

planform 2 

flat bottom 1 


3 

1 

1 

2 

1 

1 

1 


3 1 1 
3 1 1 
3 1 2 
3 1 1 
3 1 2 
2 1 1 
2 1 1 


Total CPU time for Euler design: 36 hours 


9 

3 
6 
6 
6 

4 

2 

36 


Table 3. CPU Time for Parametric Study 
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Figure 1 . F-5 A grid with H-0 grid topology used in earlier computations 



Figure 2. F-5A grid with C-0 grid topology for front block and H-O grid for rear block 


29 




a) entire crossplane grid at FS 14.02 



m 

llllllliminl 

Sill! 



M 


b) near body details at FS 14.02 
Figure 3. F-5A forebody grid details in crossflow plane 
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a) entire crossplane grid at FS 29.61 



b) near body details at FS 29.61 
Figure 4. F-5A forebody grid details in crossflow plane 
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Figure 7. Comparison of F-5A in viscid surface pressures on the leeward plane 
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Figure 8- Sign convention for forces and moments used in the present study 
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Figure 14. F-5A forebody turbulent stagnation pressure contours for a = 40° and /? = 5 C 
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Figure 15. F-5A forebody pressure vectors at FS 14.02 for a = 40° and (3=5° 
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Figure 17. F-5A vortex path along forebody for a = 40° and p = 5° ■ 
(turbulent stagnation pressure contours) 
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(b) closeup at the nose 


Figure 18. Erickson chine fore body longitudinal baseline grid details 
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Figure 20. Comparison of Erickson forebody inviscid surface pressures 
for H-O and C-0 grid topologies at FS 7.19 
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Figure 21. Comparison of Erickson forebody inviscid surface pressures 
for H-0 and C-O grid topologies at FS 13.56 
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Figure 22. Erickson chine forebody surface pressures at a = 30° and (5-5° 
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Figure 23. Erickson chine forebody surface pressures at a = 30° and /3 = 10° 
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Figure 24. Erickson chine forebody surface pressures at a = 40° and /3= 10° 












- Cp outward 
+Cp inward 


a) Inviscid 


b) Turbulent 


Figure 27. Erickson forebody inviscid and turbulent pressure diagrams at FS 7. 19 
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b) FS 13.56 


Figure 29. Erickson forebody inviscid stagnation pressure contours for a = 30° and /? = 5 
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Figure 30. Erickson forebody turbulent stagnation pressure contours for a = 30° and [3=5° 
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Figure 32. Surface Pressures on a generic forebody showing effect of multigridding and 
multi sequencing 






















Figure 34. Effect of varying b/a on the directional stability characteristics 
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a) a = 20° b) a = 30° 



c) a = 40° 


Figure 35. Effect of varying b / a on side force at various angles of attack 
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Figure 40. Stagnation pressure contours for a - 40° and p = 5° at x - 27.99 in. 
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Figure 41. Effect of varying chine angle on 
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Figure 42. Effect of varying chine angle on side force at various angles of attack 
for b / a = 0.5 
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Figure 44. Effect of varying chine angle on the variation of ACp at a = 30° 
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Figure 49. Effect of unsymmetrical shape factor n on the directional stability characteristics 
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Figure 50. Planform shapes used in this study 
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Figure 51. Effect of planfonn shape variation on the directional stability characteristics 
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Figure 52. Effect of planform shape variation on the side force at a 20 
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Figure 53. Effect of planform shape variation on the side force at a = 40° 
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